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■ We derive an expression for diffractive F2 structure function which should be valid at small 
x for quasi-elastic scattering on a hadron and for quasi-elastic scattering on a large nucleus. 
This expression includes multiple rescatterings of the quark-antiquark pair produced by the 
virtual photon off the sources of color charge in a quasi-classical approximation. We find 
that there is a relation between such diffractive production and inclusive processes. In the 
former, one averages over all colors of sources before squaring the amplitude, and in the 
latter one first squares the amplitude and then averages it in the hadron's or nuclear wave 

^"■j ' function. We show that in the limit of a large virtuality of the photon Q the diffractive 

, structure function becomes linearly proportional to the gluon distribution of the hadron or 

■ nucleus, therefore proving that in this sense diffraction is a leading twist effect. 



I. INTRODUCTION 



In recent papers by Buchmuller, Hebecker and collaborators diffractive processes are considered where 
a virtual photon scatters on a hadron or a nucleus, producing two jets out of a quark and an antiquark in the 
virtual photon's wave function, and a rapidity gap between the quark-antiquark pair and the hadron (nucleus) 
with which the virtual photon scatters. In the final state the hadron (nucleus) forms a cluster of diffractively 
produced soft particles. The process is illustrated in Fig. [l]a. Buchmiiller, Hebecker and colleagues Jl|-|| 
make the provocative claim that such processes may be computed in a semi-classical description of small x 
processes where the effect of the hadron is taken into account as classical sources of color charge They 
claim that such diffractive processes may be computed by first computing the amplitude of the process, then 
averaging it over the various color orientations of the source and finally squaring the amplitude to obtain the 
cross section. That way the squaring of the amplitude and averaging in the hadron's (nuclear) wave function 
are done in the order opposite to what one usually does when calculating total (inclusive) cross section. 

In this same theory the structure functions for deep inelastic scattering are computed using the 

same distribution of classical sources. The only difference is that the amplitude is first squared and then 
averaged over color. If these prescriptions are true, then there is a subtle and deep relation between the 
structure functions computed in deep inelastic scattering and diffractive virtual photon production of jets. 

In this paper, we find that we can derive this relation for a restricted class of diffractive processes with a 
rapidity gap where the hadron encountered remains intact (see Fig. |l|b). We shall call such processes quasi- 
clastic virtual photoproduction. Although it might be possible to extend our methods to the more general 
case where the hadron fragments into soft particles, we are unable to do so with our current technique. 

We shall present two types of arguments. The first is an intuitive argument similar to that invoked in 
Glauber scattering. We shall pursue the argument in the context of a large nucleus where one has a simple 
picture of the sources of color charge arising from the nucleon valence quark distribution. We will then 
later generalize the argument to a hadron where the sources are largely gluonic and arise from a proper 
renormalization group treatment of the hadronic wavefunction. 
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We shall perform this analysis when the averaging over color sources is Gaussian in the same sense as was 
discussed in fi|-|lO|], or more generally local in transverse space. If the renormalization group analysis after 
including the QCD evolution in rapidity |§Hl0| yields the weight function for averaging over sources which 
is not local, as it well might be, then our prescription would fail to impose the condition of a diffractive 
rapidity gap on the process at hand and should be understood in the sense of quasi-classical approximation 
only. 
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FIG. 1. (A) Diffractive deep inelastic scattering with a rapidity gap. (B) The same process with the hadron 
remaining intact in the final state (quasi-elastic photoproduction). The interaction is roughly illustrated by an 
exchange of some color singlet particle (pomeron). 

We derive explicit expressions for the quasi-elastic structure function of a nucleus or a hadron. This 
structure function is analogous of the diffractive structure function of Buchmuller, Hebecker and colleagues. 
We shall refer to it as F® with the caveat that we are capable at present of describing those processes which 
are quasi-elastic for the hadron. We also derive an expression for the structure function for deep inelastic 
scattering F% , showing that F2 and F® are related by an interchange of color averaging and squaring of the 
amplitude. 

We analyze the large Q 2 limit of the expression we obtain for the diffractive structure function F 2 D . The 
result is very interesting. The diffractive structure function of the hadron or nucleus turns out to be linearly 
proportional to the gluon distribution of the hadron or nucleus: F®(xgj, Q 2 ) ~ xsjG(xBj,Q 2 )- That way 
we prove that due to the multiple rescattering effects diffraction becomes a leading twist effect jH]], in the 
sense that diffractive structure function is not suppressed by extra powers of Q 2 with respect to the total 
structure function. At large enough Q 2 diffractive structure function should depend on center of mass energy 
for deep inelastic scattering, or, equivalently, on Bjorken XBj, in the same way as the structure function F2. 
This statement is confirmed by recent HERA data fj^] . 

The outline of the paper is the following: In Sect. II we construct the light cone wavefunction of a virtual 
photon. This will be useful as we shall derive the scattering amplitude by convoluting this wavefunction 
with the propagator of the quark-antiquark pair through the hadron which includes multiple scatterings of 
the qq pair on the gluons (classical sources) which constitute the hadron. 

In Sect. Ill we analyze diffractive scattering from a large nucleus and derive an expression for F%. This 
example illustrates all the essential points of the more general analysis where a hadron is treated as a sum 
of sources arising from a renormalization group treatment of the hadron wavefunction. We show how in the 
large Q 2 limit due to the multiple rescattering effects diffractive structure function behaves like a leading 
twist expression. 

In Sect. IV we generalize scattering off a nucleus to a hadron. We argue that if the weight function which 
is used to color average the sources is Gaussian the result of the previous analysis are essentially reproduced. 
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If the averaging procedure is local (but not Gaussian) we shall see how the result above generalizes. We 
have not been able to generalize to a non-local (in transverse coordinates) averaging procedure. 

In Sect. V wc summarize our results, and compare with the case of deep inelastic scattering by deriving 
expressions for the total inclusive cross-section and Fi in the quasi-classical limit. 



II. LIGHT CONE WAVE FUNCTION OF A VIRTUAL PHOTON 



Here we are going to calculate the wave function of quark-antiquark fluctuations of a virtual photon. The 
diagram is shown in Fig. |2|. Without any loss of generality we can work in a frame where the transverse 

momentum of the virtual photon is zero, so that the momentum of the photon is given by = ^q+, — 0^ 

with Q 2 the virtuality of the photon. The light cone momentum of the virtual photon q + is very large, so 
that in the spirit of the eikonal approximation we will neglect all the inverse powers of q+ in the calculations 
below. However we will keep all the inverse powers of Q 2 , therefore resumming all the "higher twist" terms. 




FIG. 2. Light cone wave function of a virtual photon. 



Using the light cone perturbation theory with the Feynman rules from |l3[ ] one can write down the following 
expressions for the value of the amplitude in Fig. ||: 

= e f WT 0^^f M-k)e^ .,v A k) (1) 

with r and r' being quark and antiquark helicities correspondingly, r = ±1, A is the photon's polarization, 
A = ±1, nif is the quark's mass and e/ is its electric charge (/ denotes quark's flavor). We defined 
z = k + /q + as the fraction of the photon's light cone momentum carried by the quark. We will calculate 
separately the cases of a transverse and longitudinal polarization of the photon. We start with transverse 
polarization ej, = (0, 0,e A ), with e A = (1 + Xi)/^/2. After plugging in the explicit expressions for u r and 
v r i a straightforward calculation yields 

(k, z) = e / k2 + m 2 + 1 Q 2 z(1 _ z) (*«-'k ' e A [Kl ~ 2z) + A] + rS r ,- r ,m f (l + r\)) . (2) 

For the reasons which will become apparent later we want to obtain the virtual photon's wave function in 
transverse coordinate space. Performing a Fourier transform of Eq. (|^) 

^ r ,(x,z) = J ^ e -^ r ,(k,z) (3) 

we obtain 
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^(x, z) = -± (6 rr ,ie x • V x [r(l - 2z) + A] + r5 r ^ r ,m f {l + r\)) K Q {xa), (4) 



b rr i (x, z) — ^ 
where a 2 = Q 2 z(l — z) + mj and a; — |x|. Defining 



<i> T (x,y,z) = ^ £ ^(x,z)^ r T A (y,z) (5) 



2 

X,r,r' , / 



we derive 

- — - rv. / _ y . v _ _ _ \ 

(6) 



$ T (x,y,z) = 27V C V ( a *?-y.K 1 (xa)K 1 {ya)[* i + (1 - z) 2 ] + m 2 X (:m)^oM) 

^ 7T \ xy 1 J 



When x = y the object in Eq. (g) becomes just a square of the transverse virtual photon's wave function. 

To calculate the longitudinal contribution to the wave function we write down an expression for the 
longitudinal polarization by requiring that e 2 = 1 and e ■ q = 0. We end up with e£ ^ = (Jq~, ^p^)- 
Similarly to the above one obtains 



= e / u2 | _ 2 | ^ 2Q^z(1 - z). (7) 

Fourier transformation of Eq. (Q) yields 



k 2 + m 2 f + Q 2 z(l - z)' 



^(x,z) = ^f8 rr .z{l - z)2QK (xa), (8) 

Z7T 



which results in 



^ L (^y,z) = 2N c Y,^ ML Wz\l - zfK Q (xa)K (ya)] . (9) 
/ ?r 

Eqs. (^) and (|^) provide us with the square of the virtual photon's wave function: 

$(x,y,z) = $ T (x,y,z) + $L(x,y,z). (10) 



III. DIFFRACTIVE STRUCTURE FUNCTION 



In the calculation of the quark-antiquark propagator through the nucleus we will make use of the quasi- 
classical approximation |]l5| , |l6f , which restricts the interactions of the qq pair with the nucleons to nothing 
more than a two gluon exchange with the nucleon in covariant gauge. Light cone (A_ = 0) gauge analysis 
of the process is a bit more saddle, but leads to the same formulation of the classical limit f| . 




oo ••• o 

FIG. 3. Diffractive quark-antiquark production in DIS. 
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As was stated above, in order to calculate diffractive (or elastic) production of the quark-antiquark pair 
one has to average the production amplitude in the nuclear wave function. After one does that all the one 
gluon exchange contributions would cancel resulting in two gluon exchange (virtual) interactions being the 
only possibility. The amplitude is depicted in Fig. ||. There the gluon lines can hook in all possible ways 
to the quark lines (see Fig. || ). Similarly to the wave functions, the propagator is calculated in the eikonal 
approximation, i.e., the qq pair is assumed to have a very large light cone momentum, inverse powers of 
which in the amplitude will be neglected. 

In order to calculate the quark-antiquark propagator we have to start at the lowest order case of one 
nucleon by summing the diagrams shown in Fig. ^. The vertical dashed line there indicates the final state 
of the qq pair. Because the color averaging is done in the amplitude and the complex conjugate amplitude 
independently, the propagators in the amplitude and the complex conjugate amplitude effectively factorize. 
That way we can do the summation just in the amplitude and the answer for the complex conjugate amplitude 
will be the same. That is why in Fig. we show the lowest order interactions only in the amplitude. 



=C3 C3 =(□)= 

ABC 

FIG. 4. Calculation of the quark-antiquark propagator. 

To evaluate diagrams in Fig. ^, similarly to J^,[^] we start by defining the normalized quark-nucleon 
scattering cross-section 

where 1 is the momentum transfer. The Fourier transformation of V(l) is defined by 

F(x) = / d 2 le~ U x V(l). (12) 



Then, as one can easily see, the diagrams in Fig. [| give the following factors: A and B give — (1/2)^(0) 
each, and C gives V(x), where x is the transverse separation of the pair. Defining 

v(x) = ±[l-V(x)] (13) 

and adding all the contributions of graphs in Fig. || we obtain 

A + B + C = -[V(0) - F(x)] = -ix 2 5(x). (14) 

Now to obtain a contribution to the propagator we have to multiply this result by pa = 1/A, to take into 
account the density of nucleons and quark-nucleon cross-section, and by the length of the path along the 
z direction of the quark-antiquark pair in the nuclear matter L. A is the mean free path of the pair in 
the nucleus. One can easily see that if we start including more and more nucleons, interacting with the qq 
pair, the answer for the propagator will be just an exponentiation of the lowest order result. Since in light 
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cone perturbation theory one has to also include the diagram without interactions, the final answer for the 
propagator of the pair in the nuclear matter will be 

6XP ( 4X~ J' (15) 

One should note here that our convention is different from [ fl5[ , since a in our case is the quark-nucleon 
cross-section, same as in p6[ . 

An important fact to note is that in the two gluon exchange, no net transverse momentum is carried to 
the source. The gluon coming in has the same transverse momentum as the gluon going out. Additionally, 
the longitudinal momentum transfer is very small, of order e~ Ay where Ay is the rapidity gap between the 
source and the quark-antiquark pair. This rapidity gap will in general be large. Therefore the source barely 
changes its momentum in the interaction. 

In the language of the Fock space wavefunctions one can see that the Fock space wavefunction of the 
target nucleus is essentially the same after the interaction as it was before. Therefore with probability close 
to one, the hadronic wavefunction is unchanged. Therefore the procedure of averaging over the sources 
before squaring the amplitude projects out the quasi-elastic piece of the interaction. In the final state there 
is a quark-antiquark jet, a large rapidity gap, and finally a nucleus with its momentum largely unchanged. 

Of course we expect that a real nucleus having small binding energies will fragment. Presumably quasi- 
elastic means here only no inelastic particle production in the central rapidity region. On the other hand for 
a real hadron, where we imagine the sources as gluons arising from a proper renormalization group treatment 
of the hadron wavefunction, quasi-elastic means no fragmentation of the hadron. 

This leaves open the question of how one treats truly diffractive processes where one allows for some 
fragmentation of the hadron into pions which have momentum close to the hadron's initial momentum. 
Presumably one can find a relationship to the quasi-elastic calculation, but we have not been able to do so. 

In order to write down the final answer for the qq production cross-section it is convenient to work in 
the rest frame of the nucleus. Analogous to |l5| ] we assume that the virtual photon (or the quark-antiquark 
pair) reaches the nucleus at light cone time r = X— = 0. Then, if we define t± as the time when the virtual 
photon splits into the quark-antiquark pair in the amplitude and t 2 as the time when the virtual photon 
splits into the quark-antiquark pair in the complex conjugate amplitude, one should consider four cases: 

(a) n > 0, t 2 > (16a) 
(6) n < 0, t 2 > (16b) 

(c) Tl > 0, t 2 < (16c) 

(d) ti < 0, r 2 < (16d) 
The corresponding cross-sections will be 



da^ 1 



jM**y.-~™«* r ,.) (17.) 



d 2 kdz 7r(27r) 3 



d 2 kdz 7r(27r) ; 



J d 2 bd 2 xd 2 ye- lW ^<P(*,y,z) exp f - X *? W 5 - j (17b) 



1 f dHMye-^-^^y,z) exp f j^VgZg ^ (1?c) 



dPkdz 7r(27r) 3 

d(j{d) 1 f J2u ,2 ,2 -<k.fx-y)*/ x ( x 2 f)(x)Vi? 2 -6 2 \ ( y 2 i(yWK 2 -b 2 ^ 
- 1 d 2 bd xd 2 ye 4K lx y > $(x, y, z) exp — —— exp 1 



d?kdz 7r(27r)3 J a v ' *\ 2A / V 2A J 

(17d) 
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with <£>(x, y, z) given by Eq. (flo|). Here R is the radius of the nucleus, which we assume to be spherical, and 
b is the impact parameter. Adding together the contributions given by Eqs. ( [l7j ) we find the expression for 
the quark-antiquark production cross-section in deep inelastic scattering on a nucleus: 



da 1 



cPkdz 



— [ d 2 bd 2 xd 2 ye-^-y^(x,y,z) 

7r(27r) J J 



1 — exp 



x 2 5(x)Vi? 2 - b 2 
2A 



1 — exp 



y 2 v(y)VR 2 - b 2 



2A 



Making use of the definition of the F 2 structure function 

Q 2 



(18) 



(19) 



, o {°T + OX) 

together with Eqs. (§) , (|) , © , we derive the following formula for the diffractive structure function 

F 2 = Q 2 J2 z f f d 2 bd 2 x [ dz {a 2 K 2 (xa)[z 2 + (1 - z) 2 ] + m 2 f K$(xa) + 4Q 2 z 2 (l - z) 2 K$(xa)} 

k{2-k) , J Jo 



exp 



x 2 £(x)Vi? 2 - b 2 
2A 



(20) 



where Z 2 = ot EM j 'clem for each quark's flavor. Here we should note that in the logarithmic approximation 
of small x (large Q 2 ) 



u(x) 2n 2 a s 



PXb 3 G(xb 3 ,1/x- 2 ) 



(21) 



(see Eq. (4.7) of as well as Jl9| ). Now a s is the strong coupling constant, p is the density of nuclear 
matter, normalized to give the atomic number A when integrated over the whole nucleus. XBjG(xsj , Q 2 ) is 
the gluon distribution in a nucleon, which we take at the lowest order in a s (see the next section). Using Eq. 
( pl| ) one can see now that in the large Q 2 limit the differential cross section of Eq. ( |l8| ) becomes proportional 
to the square of the gluon distribution function of the hadron or nucleus, reducing to an old and well known 
result |0|. 

However, the large Q 2 limit of is more interesting. Let us assume for simplicity of the calculations 
(only in this section) that our target is a large cylindrical nucleus with the radius R and the length 2R along 
the axis. Then the impact parameter integration in Eq.(p0|) becomes trivial giving just a pre- factor of irR 2 . 
Using the representation of the modified Bessel function -Ki(£) [J] 



£ f°° dt _,_e 



4./o t 2 



(22) 



in Eq. (|20|), assuming that the dominating values of z are small, performing the z and t integrations, we 
obtain 



dx 



D _ N C R 2 ^ 2 

67T f Ji/Q 2 \ x ) 



2 r 



1 — exp 



2N r R 2 



x 2 x B jG(x B j,l/x 2 ) 



(23) 



We thank D. Kharzeev for showing us this trick. 
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where XBjG(xsj, Q 2 ) is now the gluon distribution function of the whole nucleus. In the leading logarithmic 
approximation we assume that the gluon distribution is a slowly varying function of Q 2 . That way it will 
effectively play the role of the upper cutoff of the x integration in Eq. (|23|). Performing the integration in 
Eq. ( |23| ) we obtain the following expression for the diffractive structure function at large Q 2 : 

F?(xBj,Q 2 ) - ^ E Z / XB 3 G(x B j,Q 2 )- (24) 

A little more careful estimates would introduce a factor of 2 In 2 in front of XBjG on the right hand side of 
(p4|). One should note here that the exact coefficient can not be precisely fixed by this method, since to be 
able to use XBjG as a cutoff in Eq. (|23| ) one should reach the kinematic region of large gluon distributions. 
This would correspond to evolving x B jG(xsj , Q 2 ) with some QCD evolution. However, in the strict sense, 
our x B jG should be taken only at the lowest order in a s [see Eq. below]. Nevertheless our qualitative 
conclusions should not be altered by introducing evolution in the gluon distributions. 

That way we see that the diffractive structure function in the large Q 2 limit is linear in XBjG with no 
Q 2 suppression. Due to the multiple rescattering effects, the structure function becomes a leading twist 
expression, i.e., it is not suppressed by powers of Q 2 . By rcsumming higher twist terms we obtained an 
answer which is effectively leading twist. This result agrees with the previous predictions of pl] , pl| . Eq. 
( p4| ) shows that the diffractive structure function depends on energy in the same way as the structure function 
F 2 at large Q 2 . This result agrees with the recent ZEUS data . 



IV. GENERALIZING TO HADRONS 



To generalize the considerations above to the case of a hadron, we first recognize that the renormalization 
group treatment proposed in |8|-|l0|| replaces the gluon distribution of a hadron by a set of sources which 
include both valence quarks and sea gluons. If we go to small enough x, the gluon density per unit area is 
large enough so that one can use a weak coupling analysis. In this respect, the situation is entirely parallel 
to that for a nucleus. 

The essential difference is that the function which we use to average over sources changes. If we use a 
Gaussian distribution of sources, 



J§>) (25) 

then it is an easy exercise to show that the result of the last section are reproduced. To do this we use the 
fact that the factor A(y) is the charge squared per unit area [as defined in [^|] with y ~ lnir + the space-time 
rapidity, and the DGLAP equations |^2| . As was discussed in some detail by McLerran and Venugopalan in 
the propagator of a single quark through the nucleus can be represented as a path ordered exponential 
(Wilson line) 

Pexp J dx' + A^(x' + )^j , (26) 

where the field A- is the gluon field taken in covariant gauge. Direct calculation of the path ordered integral 
done in |^| yields for the propagator of the quark-antiquark pair through the hadron [see formula (119) in 

§] 

7 (x) = exp [a\C f x 2 ln(x 2 A 2 QCD ) A A(y)] , (27) 

where the no interaction contribution was included in the propagator in Eq. (|27j). A is the atomic number 
of the nucleus. Using the fact that 



AA(y) = — £ / dx' G(x', 1/x 2 ) (28) 
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with p = jtt^b? nucleoli number density in the nucleus we obtain 

7 (x) = exp I a 2 n^R 2 - 6 2 x 2 ln(x 2 A 2 D ) p I dx' G(x', 1/x 2 ) ] , (29) 



/(x) = exp (a 2 iry/R 2 -b 2 x 2 \n(x. 2 A 2 QCD ) p f 



To establish the correspondence between the result of Eq. (|29j) and the qq propagators from Sect. Ill one 
should make use of the double logarithmic limit of the DGLAP equation 

<} x Bo G{xb 3i Q 2 ) = — f 1 dx'G(x',Q 2 ). (30) 



dHQ 2 /A 2 QCD ) 



In the quasi-classical calculation of the previous section the gluon distribution function were taken only at 
the two gluon level. Noting that in that two gluon exchange approximation the gluon distribution is given 

by 



X B jG(x B j,Q 2 ) = — -hi-rj- — (31) 

n A QCD 



one can easily see that Eq. (p9|) reduces to 



/ mr 2 \/R 2 - b 2 \ 
7 (x) = exp ^ x 2 px B] G(x B] , 1/x 2 ) , (32) 




which matches the exponents employed in derivation of Eqs. (18) and (|20j) if one substitutes Eq. 
them. That way we showed that the result derived for the structure function of a nucleus [Eq. 
probably, also valid for a hadron. 

There are two assumptions here. First is that one can use the DGLAP equations and the second is that 
the distribution of sources is local and Gaussian. These are true when one is at high enough momentum 
transfer so that the density of gluon dN/ d 2 prdy per unit area is small. This distribution at large pr goes 
as 1/p 2 ". As the transverse momentum is lowered therefore the density becomes large and our assumptions 
break down. 

The greatest complication in the small pt region (p|, less than or of the order of the color charge density of 
gluon per unit area) is that the Gaussian distribution changes to a presumably non-Gaussian and non-local 
(in transverse space) form. In the exchanges of gluons between the produced jet of quark-antiquark pairs, 
some transverse momentum can be absorbed by the source. For local Gaussian distributions 

(p(x0...p(x w )> (33) 

factorizes into a product of terms such as 

(p( Xl Mx 2 )> ~ «< 2 >(xi-x a ). (34) 

The locality of this product and its translational invariance are what guarantee that the source absorbs no 
momentum. 

If there is non-locality in the Gaussian averaging, then one couples together sources at different transverse 
coordinates. There will be overall no transverse momenta transferred to the sum of sources, but there may be 
momenta transferred to the individual sources at their respective transverse coordinates. (We have a picture 
in mind where we have coarse grained in transverse coordinate. Each coordinate cell contains a separate 
source.) If this is the case, then the Fock space wavefunction of the hadron has been modified somewhat 
from its original form as in the amplitude one has redistributed the transverse momentum among its various 
constituents. 

The problem has therefore become immensely more complicated. One should note however that this can 
only affect the soft transverse momentum part of the hadron wavefunction, so that perhaps this does not 
spoil the wavefunction so much. It may also turn out that the non-localities in the renormalization group 
improved distribution function for sources is small. 

As a practical matter, the transverse momentum scale of interest for the non-linearities is presumably the 
virtuality of the photon Q 2 . If one is at large Q 2 , then non-linearities and non- locality is not important. In 
this case, the exponential factors in Eq. (^0j) linearize, and one is back to the simple two gluon exchange model 
for diffractive processes. Such a model has been derived in the literature [ p0[ from various considerations 
different than those above. 
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V. TOTAL INCLUSIVE CROSS SECTION 



In the calculation outlined above we have derived the expressions for diffractive (elastic) quark-antiquark 
pair production and structure function in DIS on a nucleus in the quasi-classical approximation [Eqs. 
(|l8|),(20|)1. Our results include all higher twist effects at the classical level. As was noted above the classi- 
cal limit in our case is understood as modeling the interactions by no more than two gluons per nucleon. 
Quantitatively that results in a^A 1 ^ being the expansion parameter. That way our calculation could be 
understood as resumming all powers of this parameter. 

We should note here that similar results were obtained by Buchmuller, McDermott and Hebecker in 
However, our expressions in Eqs. jl^),(pb]) give a more explicit power of the exponential in the propagator 
of the pair through the nucleus. 

We should also note that using the techniques outlined above and used previously in fl5|| one can also 
calculate the total (inelastic) qq production cross-section. The qq propagator calculation is a little more 
difficult. In the framework of the quasi-classical approach one has to also include one gluon exchanges, 
which leave the nucleon in the color octet state, therefore breaking it apart. Without giving any details we 
state the results |l8[ . The quark-antiquark total production cross-section is 



da 
cPkdz 



. 1 .„ I ' d 2 bd 2 xd 2 ye-^ { *- y) <S>{*,y,z) 
Tr(2Tr) i J 



1 — exp 



2A 



b 2 



exp 



2A 



exp 



(x-y) 2 £(x- y)Vfl2_fr2 
2A 



(35) 



where the transverse momentum of one of the quarks is fixed to be k and the momentum of the other is 
integrated over. This formula reduces to a similar formula derived in [Eq. (32)] only in the logarithmic 
limit, when we can neglect the x dependence in u(x). The nuclear structure function in this quasi-classical 
limit is 



F 2 =Q 2 



2N C 

7r(27r) ; 



E 



7 2 



d 2 bd 2 x 



[ dz {a 2 K 2 (xa)[z 2 + (1 - zf] + m 2 f K 2 (xa) + 4Q 2 z 2 (l - z) 2 K 2 (xa)} 
Jo 



exp 



x 2 £(x)Vi? 2 - b 2 
2A 



(36) 



which agrees with the previously derived expressions [see P|l4|] and references therein] . In the limit of large 
Q 2 Eq. (pq) provides us with the usual leading twist result 



dF 2 



dln(Q 2 /Q 2 ) 3tt 



= Z 2 f x B] G{x B] ,Q 2 ). 



(37) 



/ 



That way we have established a procedure of inclusion of higher twist terms in the nuclear or hadronic 
structure functions and cross sections in the quasi-classical approximation. In both cases of diffractive and 
total inclusive cross sections the calculation is simple and follows the same prescription. One has to convolute 
the wave functions of a virtual photon with the propagator of the quark-antiquark pair through the nucleus. 
For the case of total inclusive cross section one should first square the propagator and then average it in 
the nuclear wave function. In the case of diffractive process one has to first average the propagator in the 
nuclear wave function and then square it. 
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The relationship between the expressions ( |20|) and ( pq ) can be shown in a different way |23| . Let us rewrite 
the total and elastic cross sections for the deep inelastic scattering on a hadron or a nucleus in terms of the 
S'-matrix at a given impact parameter of the collision, S(b), as |23|] 

atot = 2 / d 2 b[l - S{b)} (38a) 



a el = J d 2 b [l~S(b)} 2 . (38b) 

Now one can see that the diffractive (quasi-elastic) cross sections of a qua rk-a ntiquark pair interacting with 
the nucleus, which was used in Eq. (EQ), corresponds to t he formula (|38b|), where as the total inclusive 



cross section used in formula ( |36|) corresponds to Eq. (38a). Eqs. (38a) and (38b) show how easily the 



quasi-elastic and and total cross sections could be obtained from one another. 
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